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Abstract 



The threshold for creep in the Bean critical state is investigated. We perturb the 
Bean state by an energy Ae. We find that no matter how small Ae is it will always 
be able to induce creep somewhere on the Bean profile. This finding has important 
consequences for the interpretation of low temperature creep phenomena in terms of 
quantum creep. 
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1 Introduction 

At non-zero temperature flux lines can move due to thermal activation over energy barriers. [1, 
2] The creep velocity v of the flux lines will accordingly be proportional to an Arrhenius 
activation factor 



where U is the energy barrier, ks denotes Boltzmann's constant, and T is the temperature. If 
all barriers U are larger than some minimum barrier Um the creep will vanish exponentially 
for temperatures smaller than Um/kB- It is in general tacitly assumed that Um is larger 
than zero. Experiments which find that the creep rate does not extrapolate to zero as 
the temperature approaches zero have been interpreted as indicating quantum tunneling of 
macroscopic vortices. [3] However, if the distribution of energy barriers Dili) has no gap, 
i.e. Ura = 0, detectable creep may be thermally induced at any non-zero temperature. 

We demonstrate by use of computer simulations that no minimum barrier exists for two 
different implementations of models of the Bean state in one dimension. In the following 
section we introduce the models. We explain how we establish the Bean critical state and 
how we measure the energy barriers. We present our results and end the paper with a 
discussion of the implications of our model studies. 



V ~ exp{ — [//A;Br} 
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2 Model 



Our system consists of particles interacting via a repulsive pair potential. The particles are 
called vortices and we denote the pair potential by Uyy{r)^ where r is the separation between 
the pair. In addition to their mutual interaction the vortices interact with a set of randomly 
positioned attractive pinning wells with the pinning potential denoted by Uyp{r)^ where, in 
this case, r is the distance of the vortex from the position of the pin. 

We consider two different sets of interactions. The first is when the potentials take a 
Gaussian functional form. The vortex-vortex repulsive interaction is 

UZ{r)=Ayexp{-{r/Ryf) (2) 

and the vortex-pin attractive interaction is 

U^,ir) = -A,eM-ir/R,r). (3) 

In the second, the potentials are given by truncated parabolas. Again, the vortex-vortex 
interaction is 

.rP( . _ f Ay{\r\ - Ryf if \r\ < Ry 

' 1 if Irl > Ry ^ ' 



whilst the vortex-pin interaction is 



^^■P^"^' - \ if |r| > R, 

The positions of the vortices are restricted to the interval [— L\. The positions of the 
pinning centres rf are uniformly randomly distributed throughout the pinned region (PR) 
of the system, defined on the two intervals [— —a] and leaving an unpinned region 

(UR) on the interval [—a, a]. Vortices are entered into the UR two at a time at positions 
r = ibA, with A < a. The vortices leave the system when their positions become larger 
than Z, or more precisely when \r\ > L. The PR corresponds to the actual interior of the 
superconductor, whilst the mutual repulsion of the vortices in the UR, [—a, a], mimics the 
effect of the magnetic pressure set up by the external field which makes the vortices enter 
a real superconductor. Working with a symmetric interval allows one to avoid the problem 
of defining some additional boundary force in order to push vortices into the system [4]. 
The number of vortices Nyit) obviously changes with time t, whereas the number of pinning 
centres Np is fixed for a specific realisation of the random potential. 

At a given instant the potential energy of the system is given by 

Epot = - E Uyy{\r. - r,\) + Y.U.A\r^ " r]\). (6) 
i^j 8=1 i=i 

The model is driven in the following way. Vortices move according to the over-damped 
equation of motion 

dr, dEpot 
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Starting form a vortex free system, vortices are successively injected into the UR, [—a, a]. 
After each addition Eq. 7 is iterated until all of the vortices come to rest. We say this 
situation has been reached when the largest force on any one vortex is < 10-^M„/i?p. Only 
when all motion has ceased are new vortices injected [5]. 

We continue this procedure until vortices starts to leave the system at the two edges 
r = ±L. A density gradient is set up because of the competition between the repulsive 
vortex-vortex interaction and the attractive vortex-pin interaction, with the vortices being 
most dense at the surface of the superconductor r = iba. The stability of the established 
vortex profile is probed by the following method. 

Each of the vortices are visited successively. We displace the considered vortex, with 
position To, by an amount Ax, i.e. Tq — > Tq + Ax. We then relax the system according to 
Eq. 7 and then compare the vortex configurations before and after the perturbation. There 
are two possible outcomes: 

1) The perturbation Ax does not result in any avalanche and the perturbed vortex simply 
returns to its original position, Tq. 

2) The perturbation does lead to an avalanche and a general re- arrangement of the vortex 
positions has occurred. We then measure the total displacement of the vortex profile 
given hy 6 = Yln=i{''^i ~ '"n)' where (r^) denotes the initial (final) position of vortex 
number n. 

When we displace the vortices, we only consider perturbations in the range < |Aa;| < 
1.1 Rp and the vortices are always displaced in the direction of decreasing vortex density. 
Displacements towards the high density region always result in the first outcome. Eor some 
vortices in the system, perturbations in this range only ever produce the first outcome. 
Eor the remaining vortices, a sharp threshold, Axc, exists. Eor perturbations |Aa;| < Axc, 
outcome number 1 applies and there is no instability in the system. Eor |Aa;| > Axc, the 
system becomes unstable and vortex re- arrangement occurs. It is also worth noting that for 
all perturbations Axc < |Aa;| < l.li?^, the same re- arrangement, 8, of the vortices occurs. 
An increase in the perturbation does not give an increase in the vortex re- arrangement, 8. 

Our main interest is in the properties of the performed perturbation. Eor a vortex at 
position To, we calculate the threshold perturbation Axc{ro). If no instability is induced 
for |Aa;| < l.li?^ we put Axc{ro) = l.lRp. The increase in the energy of the system 
corresponding to the perturbation Axc{ro) is denoted by Ae. As we shall see in the next 
section there is no gap separating the measured values of Ae from zero. 



3 Results 

We present our main results in Eig. 1. The most important point to notice is that the 
distribution of activation energies measured in the simulations has support all the way 
down to zero. The detailed shape of the distribution is model dependent, which is seen by 
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comparing Fig. la and Fig. lb. The larger support at low energies ol the parabolic model is 
probably an effect of the sharp cut-off introduced in this model. The reason being that the 
forces derived from these potentials depends discontinuously upon distance. Nevertheless, 
both models exhibit susceptibility to perturbations of vanishing energy. It is also interesting 
to note that barrier distributions with seemingly diverging support at small energies have 
been measured experimentally on epitaxial films. [6] 

The inserts in Fig. 1 show the spatial dependence of the size of the activation energies. 
We have plotted Ae as function of the position Tq of the vortex under consideration. One 
sees that the majority of avalanches are released due to perturbations at the foot, |ro| ~ 
of the pile in the case of the Gaussian model. The spread in the measured values of Ae is 
also much larger for vortices placed in the bottom of the pile than for vortices placed at the 
top. Again the model with the cut-off is different. The response of the vortex profile in this 
model is much more spatially homogeneous with a band of Ae values for all positions in the 
pile. 

Fig. 2 further illuminates the spatial dependence of the stability of the pile in the 
Gaussian model. The plot contains a measure of the average displacement needed in order 
to induce a perturbation as function of position in the pile. A word of caution is needed. 
As we mentioned above Axc{ro) is set equal to l.li?^ if no rearrangement in the pile was 
induced as a result of displacing a vortex at position Tq a distance up to l.li?^. When the 
graph in Fig. 2 assumes the value l.li?^ this simply means that no disturbance of the profile 
was ever induced by displacing a particle at that position. In reality the average value of 
Axc at such a position might of course be larger than l.li?^. The dip in {Axc{ro)) at values 
of To close to L shows that the pile is most sensible to perturbations in this region. 

The distribution of induced creep is shown in Fig. 3 for different levels of the activation 
energy. As a measure of the creep caused by a perturbation we measured (as described 
above) the total amount of displacement 6 of the vortex positions. The distribution P{S) 
depends not very strongly on Ae. 



4 Discussion 

Our main results of our model study follows: 

(A) The stability of the vortex profile in the Bean state is strongly position dependent. 
The low density region is much more unstable than the high density region 

(B) No matter how small the energy of a perturbation is it will be able to produce motion 
in the profile. 

Since our numerical models are qualitatively similar to the physical Bean state in a real 
superconductor it is worthwhile exploring the relevance of our findings. 

Let us first describe an experiment which will be able to assess the relevance of the 
spatial dependence of the stability. Apply a temperature gradient to the Bean profile in the 
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following way. Imagine we have established the Bean state by increasing the external field 
up to a certain value. Let the temperature in the centre of the bulk of the sample be Tj, and 
the temperature at the sample surface be T^. We fix the temperatures at the values = To 
and Tfc = To + (5, where 6 <^ Tq and To is some appropriate low temperature. In this case 
we measure the creep rate Si. Next we fix the temperature in the bulk at Tj, = Tq and the 
temperature at the surface at = Tq + 6. We measure creep rate 5*2. If result (A) applies 
we will find Si ^ 82- This is because the most thermally active region coincides with the 
most unstable region when the temperature increases in the direction of decreasing vortex 
density. 

The specific distribution of perturbation energies shown in Fig. 1 might not apply 
to the superconductor for several reasons. First, our models consist of particles in one 
dimension interacting through a short range pair potential and not three dimensional fiux 
lines with long range interactions. Secondly, we only consider a very special subset among 
all possible perturbations. In reality thermal perturbations will involve the displacement of 
many vortices simultaneously, whereas we simply displaced a single vortex at a time. And 
finally, our pinning centres are relatively long ranged compared with point pinning centres in 
superconductors of range the superconducting coherence length. Although these differences 
might indeed produce a barrier distribution Dili) with another functional form than the 
one in Fig. I, we do not see any reason why any of the mentioned factors should introduce 
a gap in Dili). 

Let us finish with a discussion of the possible consequences of the lack of a gap in Dili) 
for the behaviour of creep at low temperature. One can adapt a simple thermally activated 
fiux motion (TAFM) analysis as in Ref. [7]. The creep rate {S)u averaged over the barrier 
distribution is given by 



Let us assume in accordance with our finding (B) above that the behaviour at small U can 
be described by Dili) ~ U~" and that Dili) effectively vanishes (becomes exponentially 
small) for U > Umax- The low temperature behaviour of the creep rate is then given by 



Thus, if the barrier distribution for the superconducting Bean profile behaves like Dili) ~ 
IjU for small U (in qualitative agreement with Fig. lb) we find that the thermal creep 
rate approaches a constant as the temperature is lowered towards zero. This is a behaviour 
often observed experimentally [8, 3], though traditionally ascribed to macroscopic quantum 
tunneling of vortices [8, 2]. 

We can not claim that the behaviour of our simple one dimensional models proves that 
creep observed in the limit of zero temperature is thermally activated. Additional simu- 
lations in higher dimensions and of more realistic models are needed. Nevertheless, we do 
believe that our study shows that the perturbation spectrum of the Bean critical state is a 
subtle quantity. It is difficult to measure the energy barriers directly experimentally. One 
often has to rely on deconvolution procedures involving various assumptions concerning the 
nature of the time dependence of the magnetization. [8] This might explain why the barrier 






(9) 
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distributions deduced from creep experiments sometimes exhibit a lower energy cut ofF[8, 7] 
and sometimes not [6]. Experiments designed to measure the nature of the barriers in a 
direct way are clearly of great interest. 
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Fig. 1a 
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Fig. 1b 
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Figure 1: A) The distribution of activation energies, P(Ae), for the Gaussian interaction 
modeL The model parameters are Ay = Ry = Rp = 1 and Ap = 0.5. The inset shows 
the spatial variation of the activation energies, where Tq ~ is the surface of the sample 
and To ~ iv = 120 is in the sample bulk. B) The same as in A) except for a model with 
cut-off parabolic potentials. The same values of the model parameters are also used. The 
distribution behaves like 1/U for small U. 
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Figure 2: The average displacement of a vortex at position Tq needed to cause an instability 
in the system, as a function of the position Tq. See the text for a more comprehensive 
explanation. 

Fig. 3 
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Figure 3: The distribution of the total displacement of the system after an instability has 
been induced. The two curves are for different ranges of the activation energies, Ae. 
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